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ABSTRACT 


A theory for buckling of three-dimensional space frames by 
finite element methods is developed. A procedure for automation of 
buckling calculations is presented. An accuracy study is made of the 
convergence of the theory by f ini ; e difference methods and applied to 
the pinned, fixed, and cantilever column problems. The semirigid 
joints are introduced by a linear spring approach. Two— and three- 
dimensional numerical examples are presented that show the influence 
of semirigid joints on buckling loads. The procedure used to solve for 
the buckling loads in the examples is determinant plotting. The results 
for the problems studied show that practical moment resisting joints 
may he assumed to be rigid with a slightly nonconservative buckling 
load. Also, it is shown that practical moment free connections have a 
substantial amount of residual strength for buckling. 
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CHAPTER X 


INTRODUCTION 

The buckling of planar frames has been the subject o* an 

( t_ p) 

extensive amount of experimental and theoretical study. * ' A 

survey of these investigations^) has led to simple formulas from 

which buckling loads can be calculated. Most buckling analyses 

assume that the joints are rigid. For real structures, however, the 

joints are not 'igid but have some flexibility due to the method of 

connection (bolted or welded). Such practical joints can be termed 

"semirigid. ■' If it is assumed that the joint is rigid and is 

actual semirigid, the buckling load will be overstated. A joint 

assumed to be pinned and is actually semirigid, the buckling load 

may be understated. A safer and more economical design would result 

if the effect of the semirigidity of the connections were included 

in analyzing a frame. Some work on the influence of semirigidity 

for frames has been accomplished for planar frames. p Q r 

instance, C-oering^ 1 ) used a generalized slope deflection method 

with fictitious members to represent the semirigid joints. Monforton 

( 2 ) 

and Wu modified the stiffness matrices for the members by using a 
correction matrix. Also, the methods used for planar studies (rigid 
or semirigid) do not appear to be easily extendable to three- 
dimensional buckling and another approach seems to be required. 

An approximate method which has merit for handling three- 
dimensional buckling including the influence of joint flexibility is 


1 


2 


the finite element method of structural analysis. This method is a 
Rayleigh-Ritz approximation as applied to the displacement shapes. 

It is well suited for calculation on a digital computer and has the 
potential of being applicable to almost unlimited cases of complicated 
structures composed of beams, plates or shell segments. Because of 
its advantages the finite element method has been chosen for use in 
this study* 

The purpose of thes thesis is to extend the finite element 
method to compute buckling loads of three-dimensional structures 
allowing for joint flexibility. A linear spring is used to 
approximate the semirigidity of the connections. Torsional stiffness 
for any given member is provided; however, torsional semirigid 
joints are not considered. Twist buckling is not included in this 
study and limits The sample problems to frames composed of cross 
sections consistant with simple beam-column behavior. The derivation 
and a study of accuracy of the three-dimensional buckling matrices 
are presented, as well as sample problems for two- and three- 
dimensional frames having both rigid and semirigid joints. The 
development assumes inextensional buckling behavior for all compression 
members. The stiffness values used for the semirigid joints were 

(5) 

based on the series of tests performed by Rathbun on 
representative joints for real structures. 



CHAPTER II 


MATRIX ANALYSIS BY FINITE ELEMENT METHODS 

2.1 Derivation of Three-Dimensional Stiffness and Geometric 
Matrices 

In the finite element method it is customary to obtain a 
characterization of the stiffness properties, {kJ , of each element and 
relate the end nodal displacements, <«> , to the end internal forces 

® by ( F ) - [ K ] {>} • u> 

A symbolic finite beam element with ends at points 1 and 2 is shown in 
Figure 2.1.1. The vectors /uA and (f A are the respective vectors of 

| 1 ^ Z ___ 

fi) > i*0 

Figure 2.1.1- Finite element. 

end displacements and force resultants. Thus Equation (l), as applied 
to Figure 2.1.1, is 

f Ft 1 __ i 

(2) 


W 


where [kJ is the stiffness matrix of the finite element. For buckling 
problems there are two parts to [ k] (ref. 6) in Equation (l), the bending 


stiffness portion {k°j and the portion dependant on P ^K^J. Thus, 
for buckling. Equation (2) can be written in the form 







( 5 ) 


three-dimensional column buckling. 
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k 


The sign convention for rotations used in this thesis follows 
the right hand rule and is shown in Figure 2.1.2. Consider a 
continuous prismatic beam element with a length h and a symmetrical 
cross section subjected to a constant axial load P. The strain energy 
of the beam which takes into account bending about both axes and 
constant torque is 

0 ■ | J* [VU + ■v?** * g «*,*] w 

The terms El, GJ are the bending and torsional stiffnesses and v, w, 

$ x are the relative displacements of the beam element (Figure 2.1.2). 

The distance along the element is x and the subscript x following 
a comma denotes differentation with respect to x. If inextensional 
behavior is assumed, the bowing effect of a differential segment dx in 
the x-y and x-z planes of Figure 2.1.3(b) due to out of plane motion is 
e = dx - ^dx 2 - dv 2 - dw 2 )^ 

e “{} ■ t 1 ■ (f,x * J 

Expand the total bowing term of Equations ( 5 ) in a binomial series. 

The result for the segment is 



The external potential energy of the bowing effect for the element is 


n 



e 


Sc dA 


Sc = 


p 

I 


where 



9 

y 



» 


9 

x 


Figure 2.1.2.- Sign convention. 



Figure 2.1.3.- Bowing effect. 
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Then, if the axial force and the area of the member are constant: 

” ‘I Jn f T ,x 2 h) 


The total potential energy of the element is the sum of Equations 
(4) and (7): 

v = n + u 




+ El 

y 




e 


x,x 


+ p 



+ w 

} 



dx 


( 8 ) 


The Rayleigh-Ritz method is used to approximate the behavior 
of the finite element. The displacements shapes (v,w) are assumed 
that satisfy the boundary conditions at each end. The simplest set of 
approximations that satisfy the boundary conditions - the arbitrary 
displacements and rotations at the ends of the element - is the 
following: 9^ = a^ + a-^ x 

v = b 0 + ^ x + b 2 x 2 + bj x^ 

w = Cq + c-^ x + C 2 x 2 + Cj x3 (9) 

where a^, bj_, and c^ are constants of integration. These displacement 
states correspond to that resulting for a prismatic beam undergoing 
small displacements and subjected to end loads. 

It can be seen from Equation (8) and Equations (9) that the 
bowing effect as well as the varicus energies are uncoupled from each 
other and can be handled separately. Therefore, the bending derivation 
will be restricted to two dimensions to compare with Martin’s work 
(ref. 6) and these results will be subsequently utilized to obtain the 


7 


complete three-dimensional behavior. 

Solving for the constants of Equations (9) in terms of the nodal 
displacements of the element results in the following: 


®x “ (®X2 “ ®xi^ x 


(10a) 



Let V-^ be the two- dimensional portion of the potential energy and Vg 

the remaining portion. Using the calculus of variations (incrementing 

by a virtual displacement). Equation (8), excluding v and w^ x terms, 
is 

v i+ AVi- 1 j * Tei, + sj 

Expanding and grouping together those terms linear in the virtual 
displacement 6V 1 and quadratic as 8^/^ yields 
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V 1 +AV 1 = \+ S\+ 8\ (12) 

Prom energy principles it is known that 6V = 0 for equilibrium and 
6% > 0 for stability (See ref. 7). Since the only term that is of 
interest in buckling is the second variation & 2 Vi, it will be evaluated. 
Note that for this problem the potential function is quadratic and the 
variations higher than the second vanish. Thus 

Al - 5 jC(S 6 ' T « * 7 ( + P 6 ' v} “ <«> 


The appropriate variations of Equations (lQa,) and (10b) axe substituted 
into Equation (15) • The result after expansion and integration is 
the following: 

i|lr (***& 



Equation (l4) can be written in matrix form as 


&2 Vl =|{6u) T [K] (8u) 


( 15 ) 


Separating K into two matrices, K 
result : 



matrix. 
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and K' gives the following 


6v_ 60 60 

2 x 2 « 2 



is shown along the top of the 
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6v, 69 60 5v_ 60 60 

1 *1 *1 2 *2 *; 



Equation (l6) is the standard, two-dimensional finite element 


stiffness matrix and Equation ( 17) is the two-dimensional geometric 

( 6 ) 

matrix repoi ced by Martin, but derived in a different manner. 


To extend this in the third dimension. Equation (10c) and the 
remaining portion of Equation (8) are used to obtain, through the same 
procedures as Equations (ll), (12), (13), and (l4), the following 
additional terms to the second variation (6 2 V2): 
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It is seen in Equations (l6) and (17) that there is no axtenr- 
sional contribution in the K° and K' matrices. Extensional behavior 
is neglected in the buckling behavior of the element corresponding to 
the usual inextensional assumption made when deriving the differential 
equation for column stabily. The final three-dimensional finite 
element stiffness matrix K? obtained by adding the applicable part of 
Equation (18) to Equation (l6), is the following: 


BV 1 

>? 

CD 

CO 

£ 

CD 

oO 

8 \ 

^2 

6 “2 

80 y 

y 2 

80 

*2 

12EI 




to. 

12EI 




6el 

z 

h^ 

0 

c 

0 

z 

h 2 

z 

i? 

0 

0 

0 

z 

h 2 


12EI 


6ei 



12EI 


6EI 


0 


0 - 


0 

0 

"T* 

0 

‘if 

0 

0 

0 

GJ 

r 

0 

0 

0 

0 ■ 

GJ 

‘h 

0 

0 


to: 


to 



6ei 


2EI 
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‘if 

0 

l 

h 

0 

0 

if 

0 


0 

6EI 




to. 

6EI 




2EI 

z 

h 2 

0 

0 

0 

z 

“E 

Z 

"T 

0 

0 

0 

z 

— 

12EI 




“ to. 

12EI 




6EI 

z 

"7“ 

0 

0 

0 

z 

h 2 

z 

h^ 

0 

0 

0 

z 

’IT 


12EI 


6EI 



12EI 


6EI 


0 


0 

if 

0 

0 

__Z 

h* 

0 


0 

0 

0 

GJ 

T 

0 

0 

0 

0 

GJ 

h 

0 

0 


6EI 


2EI 



6EI , 


to 


0 

■if 

0 

IT 1 

0 

0 

~f 

0 

I 

h 

0 

to. 




2EI 

6EI 




to 

T 

0 

0 

0 

z 

h 

z 

h 2 

0 

0 

0 

z 

«l 


( 19 ) 
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The final three-dimensional finite element geometric matrix 
K* is obtained by adding the rest of Equation (l8) to Equation ( 17 ): 


K = P 


6v 

5v 

60 

60 

60 

6v 

6v 

60 

60 

60 

1 

1 

*1 

y i 

Z 1 

2 

2 

X 2 

y 2 

*2 

6 




1 

6 




1 

5h 

0 

0 

0 

10 

‘ 5h 

0 

0 

0 

10 


6 


1 



6 


1 


0 

5h 

0 

" 10 

0 

0 

“ 5h 

0 

" 10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1 


2h 



1 


h 


0 

" 10 

0 

15 

0 

0 

10 

0 

" 30 

0 

1 




2h 

1 




h 

10 

0 

0 

0 

15 

" 10 

0 

0 

0 

" 30 

6 




1 

6 




1 

' 5h 

0 

0 

0 

" 10 

5h 

0 

0 

0 

" 10 


6 


1 



6 


t 

J- 


0 

" 5h 

0 

10 

0 

0 

5h 

0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1 


h 



1 


2h 


0 

" 10 

0 

‘ 36 

0 

0 

10 

0 

13 

0 

1 




h 

1 




2h 

lo 

0 

0 

0 

" 5C 

- 10 

0 

0 

0 
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( 20 ) 


The equilibrium equations describing the end forces in 
terms of any set of nodal displacements (virtual or real) is 


{f} = [K° + k'] (u) 


( 21 ) 


where 
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2.2 Introduction of Semirigid Joints 

In this study the influence of semirigidity will be studied. 
Some method of introducing joint flexibility must be provided. The 
method selected is the use of a linear spring between members at a 
joint. The strain energy of the spring is 

r 9 

U = / M d9 where M = R0 
s Jq s s 

U m f R0 dO (23) 

s J Q 


U 



Where R is t..e spring or joint stiffness. 


Some typical examples of this approach for semirigid! ty 
are shown in Figures 2.2.1 and 2.2.2. In ell cases the column angle 
is used as the base angle (datum plane) for determining the slippage 
of the Joint. A typical deformed joint is shown in Figure 2.2.1(a). 

The total strain energies for the two-dimensional springs of 
Figure 2.2.1 (a) and (b), respectively, axe as follows: 

u s - r ( 9 i - 9 2 f < 2l *> 



(24b) 


The total strain energies for the same joints hut adding the 
third-dimensional member, Figure 2.2.2(a) and (b), respectively, are 
as follows: 


I - r ( e 6 - e 7 ) 2 + £ ( e 6 - *i f 


I u s * t (°8 - \T *r( e 8- 9 io f * r ( 9 8 - * s f ^ 


For example, the matrix format of the second variation terms 
that come from Equation (2Va) is 



where R® is the spring stiffness matrix for the joint of Figure 2.2.1(a). 
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The semirigid connection may be thought of as a locally- 
weakened section between the end of a beam and the face of the column 
to which the connection is made. Reference (l) points out that only a 
small part of the rotation of a semirigid joint is due to bending of 
the angles. Most of the rotation is due to slip and extension of 
its rivets or bolts. Since this semirigid joint is quite complex 
and involves other types of movement such as local bending of the 
column, this thesis will be concerned only with the use of a linear 
rotational spring to represent the semirigid joint. 

The rotational springs of the column and beam have additional 
nodal stiffnesses,^) coluinn and (^p) beam * The rigid joint 
stiffness form can be expressed as the following: 


column + ^l^beam 

4ei 4ei 

h column h beam 


(27) 


Equation (27) must be included in Equation (26) to obtain the 
complete semirigid stiffness of the joint. The result is 



(28) 
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where[sf] is the semirigid joint stiffness matrix of Figure 2.2.1(a). 
2-3 Automated Technique of Matrix Methods 


2 . 3.1 Procedure for stress analysis . The process for 
generating a computer program of any structure that is composed of 
many finite elements is to first pick a set of coordinates convenient 
for a typical element. The generalized element displacements axe <?> 
and forces axe <?> . The energy in terms of element coordinates is 


The displacements 
sponding to ends 


u - | 0) T 00 (?) 

and stiffness matrix [VJ 
i and j 


(29) 


axe partitioned corre— 






(30) 


(31) 


In the process of connecting elements, it is found that one 
element's local coordinates are not the same for another element. 
Therefore, a set of system coordinates is chosen that is convenient 
for a system of elements and the local coordinate points are numbered 
(points 1, 2, 3; • • • ) . 

A systematic numbering process for the node points and members 
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is chosen. The element properties (E, A, I, L) are given. The stiffness 
for each element is calculated in local coordinates where i and J 
refer to the end points of each element. 

If the system coordinates are called , the transformation 
from an element's coordinates to a system's coordinates is accomplished 

/ o\ 

by a transformation matrix, |^. ' ' That is 

(H#> (32) 


The stiffness of the element is transformed to system coordinates by 


Equation (32) . 


owe#] 


Consider severe elements that are connected. The next step 
is to generate the master stiffness matrix by summing all member 

stiffnesses. It can be shown that the elements of M can be 
obtained as follows: r n r-» r ~\ 


M . ■ l[\ 


All 

members 

For all displacements that are constrained(boundary conditions) , 
it can also be shown that this can be handled by striking out rows and 
col umn s corresponding to the displacement constraints, resulting in the 
reduced master stiffness matrix 

Finally, the system applied external forces {q} are identified 


and the equation 




is solved. 


2.3.2 Buckling procedure. For the buckling 
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problem, the steps are essentially the same as for the preceding 
section except that the stiffness matrix K for an element is composed 
of two parts, K°, the stiffness of the structure and K', which 
contains P — the eigenvalue portion of the problem. The remaining 
steps are the same with the external concentrated axial loads at 
the joints equal to the internal forces of the column elements in 
the K' matrix. Since the external lateral forces and moments at the 
junction of the beam elements are equal to zero in all cases studied 
in this thesis, the equilibrium equation for inextensional behavior is 



( 56 ) 


A nontrivial solution to equation (36) requires that the determinant 


K° + K' 
R 


R I 


= 0 . 


( 37 ) 


A similar condition results for stability based on the second 
variation of the potential energy of the structure. This is then 
solved by some standard procedure for d, termining eigenvalues such 
as determinant plotting, matrix iteration, or other procedures. The 
method used in this study is determinant plotting. 

When the eigenvalue has been computed, the mode shape or 
eigenvector is obtained by assuming a value for one of the components 
of the eigenvector and solving for the remaining components in terms 


of it. 


2 . 3.5 Semirigid jointB . The process of allowing for joint 
flexibility adds only one additional step from that described in 
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the previous sections. After transforming the element coordinates 
into system coordinates, the next step is to add the system spring 
stiffness matrix (Rfj3 similar that in Equation (26) 



(58) 


to the transformed element matrix or matrices. For instance, if a 

semirigid joint occurs between points i and j, the resulting 


stiffness matrix is 


P ij] - [ K u + B ±j] <» 

Any combination of Sections 2.3 .1, 2.3 .2 and this section can 
be used to give values (two- or three- dimensions) for stresses, 
displacements, eigenvalues, or eigenvectors either by allowing for 
joint flexibility or by using rigid joints. 


CHAPTER III 


STUDY OF CONVERGENCE 


3.1 Introduction 

Since the finite element method is an approximate numerical 

solution to structural problems, some knowledge of its reliability In 

these applications is desirable. This section gives the resulte of 

an investigation of the convergence and accuracy of the finite element 

matrices derived in Chapter II. The results are limited to simple 

planar buckling; however, they provide a basis for extending the 

methods to more complice ted problems. There sure many ways to study 

convergence and accuracy, most of which are based on numerical 

solutions to special problems. The approach used here is that of 

classical order of error analysis procedures borrowed from finite 
(9) 

difference methods y*' that is, to expand the discrete equations in a 
Taylor series about the i point. 


o 

3.2 Order of Error Analysis for K + K' 

Consider a continuous prismatic beam element of length h 
subjected to a constant axial load P (Fig. 3*2.1). The generalized 
displacements are taken to be rotation 8 and deflection v at the 
end of the element. 


P 


J-— P 




P 



Figure 3*2.1.- Beam element. 



(b) 
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Since the out-of -plane displacements and rotations are uncoupled from 
the in-plane displacements, it allows for separation of the two and 
only in-plane behavior will be considered to evaluate buckling 
reliability. 


If two elements are Joined together at the i^ point, two 
finite element equations result from the sum of forces and moments 
at the point. They are respectively, 



To study the convergence character, the terms in these e .ons 
are expanded about the i th point by using Taylor series such as 

± ' , h 2 " ± 

v i± i = v ± i liv ± + 27 v 1 i • • • , 


This leads to 




2 1 1 
h Vi 
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and 



h 3 ... 

y v i 


5 

So V : 


h 7 

366(7)' 



(¥ - 3) («i 


+ h 8. 


h 4 

15 


IV 


h 6 

3^o 


VI 


h 8 .. 

55of5ST e i 


VIII^ 
1 ) 



Since 9^ Is not an Independent variable in the continuum sense it 
is useful to eliminate it from Equation (42) and (43). By solving 
Equation (43) for 9^, obtaining the successive derivatives and 
substituting these derivatives back into Equation (43), it becomes 


h ■ v i + (jho ‘ riir] + f ( h6 ) = 0 


(44) 


Substituting the derivatives of 9^ into Equation (42) leads to 


v iv . Hl- + h 4 

V i El + n 


VIII 

i P VI 

720 350EI V i 


■o(h 6 ). 


(45) 


Equations (44) and (45) show the convergence character of the 
finite element Equations (4o) and (4l) and show that they converge 
to 


and 


9 i - v i a 0 

t I 

17 n 

v i ■ TST" = 0 


(44a) 

(^5a) 


with both equations having an error of order h 
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3*3 Order of Error Analysis for K° • K* for String 

Reference (6) suggests that a simpler method he used to obtain 

K* by the superposition of a string and a beam without any interaction 
between the two. This combination then leads to the same result of 
K° previously derived but a modified K*. The K f for the string 
as applied to the planar column is 


P 

h 


5v i 

i 

0 

-1 

0 


89 , 


8v„ 


69, 


0 

0 

0 

0 


-1 

0 

1 

0 


0 

0 

0 

0 


(k6) 


By applying the same procedure as before, the two finite element 
equations at the i** point are 

^jr + S) (- V 1-1 + ^1 - T rn) (- e i-i + Vi) ■ 0 
(fr) (- Vi + Vi) + ¥ ( 9 i-i ♦ k \ + 9 i + i) ■ 0 


Expanding as before results in 


nr_ 1,2 " ^ 


.8 


„ iv h 6 „ vi h u „ vni\ 

12 v i " 5S0 V i " 560(567 V i J 


+ f^Y 2h 9' 

. V nJ\ 


h' > til h 5 y u7 VTTl 

+ r 9 i + 6 o 0 i + 7(5S6T 9 i 7 =0 


(W) 


and 
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„ ' h 5 ' " h 5 V iJ VII 

2h V i - ~ V i " 56 V i “ 2520 V i , 


' t lA TV V.6 ITT u8 


. 2EI /, .2 ' * IT IV h w . VI X h' 

+ ~h" l 6e i + h 9 i + 12 e i + 55o 9 i + 350 


e i mi - 0 


Solving Equation (49) for 9^ and substituting successive 
derivatives back into this equation leads to 

). J> 7111 

' h 4 V h v i 

9 i " v i + ibo v i + ~ w%7 ' = 0 ( 


Substituting the derivatives of 9 into Equation (48) and collecting 


terms gives 


, IV ^i 2 ^i 
i " ~W~ " 12E] 


v. rai Pv^ 


, v.2 i . J* v i "i . ^6\ „ 

h " 12EI h " 720 " 3S0ElJ )“ 0 


Equations (50) and (51) show that the string approximation 
converges to the required continuum equations but Equation (51) has 
Thus its expected accuracy is less than that given in the 
previous section. 

3.4 Application of Order of Error for K° + K' 

The results on order of error can be applied to obtain some 
further information on the accuracy of the discrete approximations 
in a particular set of problems. (10,11 ) Consider a pinned column 


with effective length, L ft , and subjected to a constant axial load, P. 
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If only the first order of error term is retained, the differential 
equation representing the finite element equation is 






VIII, 


+ C, 



= 0 


(52) 


where 


*r 


p_ 

EI 


The constants and taken from Equation (^5) have 

11 4 

values of - 611(1 respectively. Since h is small. 

Equation (52) can be solved utilizing the perturbation procedure in 

references (lo) and (ll)* Let 

\ = 7 ^ + h 4 (55) 

Where corresponds to the order of error term. 

Substituting (53) into (52) yields (retaining up to terms) 

VIII 


V* - V] + « 8 * V* - \ V] - o 


Assume 


The first part of Equation ( 54 ) becomes 


(5*0 

(55) 


*>fe' 


and this raust be zero independent of hr terras, thus 

.2 




( 56 a) 


( 56 b) 
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^ = A L 


2 


\ = - (**) £ 

The second part of Equation (5^) becomes 


■.-A® -At© •■>(<*) «©" 


Therefore 


T _ U*)‘ 

l - ” 720 


2 , Jl- 


\--w 2 * ♦£)**?£ 


Now let N equal the number of elements per half wave, this is 


The final equation form is 


P CR. = " 1 + 
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Where the negative sign says that for buckling the member must be 
in compression. This shows that the error in a buckling load calcula- 
tion is determined by the contribution of the principal error term 


e, where e = 


N 4 720* 

of elements are as follows: 


Representative values of € for various number 


20 $ 

10 $ 

1 $ 


0.91 

1,08 

1.92 


(61) 


Then, two elements per buckling half wave gives the buckling load to 
about 1 $ error. 

3.5 Application of Order of Error for K° + K' for String 

The came procedure as the previous section is used with 
o 

order of error h . The differential equation of the column effect 
approximated by a string is 


tv • ' P 

y ± " \ v i + “ 


C X v IV 
°1 i v i 


= 0 


(62) 


where 


\=\ + h \ 

The constant =- from Equation (51) and Eq. (62) becomes 




= 0 (63) 


Solving this in the same manner as before 


p s - (nut ) 4 
CR. v 7 


1 + 


1 it 


El 


( 64 ) 
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JT 1 " 

The principal error term e = -o — 

NT12 

as follows: 

€ 


for various number of elements are 


N 


20$ 

2.03 

10$ 

2.87 

1$ 

8.72 


Thus nine elements per buckling half wave are required to obtain the 


buckling load to about 1$ error. 


3*6 Numerical Comparison 

To complete the study of convergence, numerical results for 
the classical columns (pinned, fixed, and cantilever) were obtained 
for several numbers of finite element segments per column. 

The results for Figures 3*6.1, 3*6.2, and 3 . 6.3 were obtained 
by the o( Infinite element approximations (K° + K f ) and shown in 
Table 3. 6.1 (a). If it is assumed that the numerical solution for 
V approaches the exact solution according to the order of error, 

(12) 

the exact result can he extrapolated hy assuming a linear relation 
for order of error (h/L)\ Also, the perturbation solution to the 
finite element equation can he applied to simple problems. It is 
given in Chapter 3, Equation ( 60 ) and gives higher but reliable 
estimates for the results given in Table 3 . 6.1 (a). The perturbation 
results were obtained by solving Equation (60) . 

Results were also obtained for the K° + K' for string, 0(h 2 ) 
for the case of a pinned column. As an example, the data is given in 
Table 3*6.2 and plotted in Figure 3*6.4. The numerical data in 

p 

Table 3*6.2 also shows a close comparison between the 0(h ) numerical 


P C R. l 2 /ei 


Number of elements 


4 ! EXTRAP 


Exact 


Pinned 12. 0000 | 9* 9200 | 9-8847 


Fixed 


9.3760 9.8696 


4o.oooo 40.3432 39 . 7754 39.5126 39.4784 


Cant. 2.5000 2. 4855 2.4711 1 2.46761 2.4674 


( a ) Numerical calculations 


P CR. L / EI 


Number of elements 


Pinned 


Fixed 


2 

3 

9.9531 

■ 

2.4883 

2.4715 

44.8194 

40.5334 


Exact 


(4>) Perturbation solution estimates 
TABLE 3*6.1.- N0NDIMENSI0NAL BUCKLING LOADS, Q(hl FOR K' 

_ __ P CR. L 2 /EI j 


Pinned 

Number of elements 

77V 0 r» + 

col. 

2 

3 

4 


Num. 

12.0000 

10.8000 

10.3866 

9.8696 


11.8990 

10.7715 

10.3769 

9.8696 


TABLE 3.6.2.- NQNDIMENSIONAL BUCKLING I©ADS OBTAINED FROM 
"STRING" APPROXIMATION TO K 1 - ©(h) 
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results and the perturbation estimates of the numerical results. 

They also demonstrate that the string finite element does have an 
2 

order of error of h . Therefore, from a practical point of view of 

seeking the greatest accuracy from the least amount of work, it is 

obviously advantageous to use K° + K* matrices having C^h 1 ) in lieu 

of the K° + K* for the string having o(h 2 ) . 

The results also show that the C^h 1 ) approximation gives good 

4 

accuracy and has excellent convergence properties. The 0(h ) 
approximation results were high but were converging to the correct 
answer as the number of finite elements was increased. The finite 
element approximation is used in the remaining portion of the study 
to obtain calculations for frame buckling with and without joint 
flexibility. 


CHAPTER IV 


PRESENTATION OF NUMERICAL RESULTS 


4 . 1 Introduction 

The homogeneous linear equations expressed by Equation (36) have 
a trivial solution, that is all displacements equal zero, and a non- 
trivial or unique solution where the determinant of the [k° + kD 
matrix equals zero. Since the trivial solution is 
useless, the unique solution is the one of interest here. By 
holding all values in the K° and K f matrices fixed, the only 
parameter is the axial load, P. The procedure used in this paper to 
determine the buckling load is to obtain the lowest nonzero value of P 
(load) that causes the determinant to vanish. This is accomplished by- 
plotting the determinant versus the magnitude of the load. The 
numerical results were calculated using the Control Data 6600 digital 
computer at the Langley Research Center. Since the results are 
obtained for three- dimensional frame structures, several types of 
member properties are required. To give realistic results, specific 
dimensions were selected for the representative structure. The columns 
are all taken as 8M34.3 steel members and the beams as standard AISC 
0WF31 sections. Table 4.2.1 shows the properties of these members. 

4.2 Rigid Joints 

Several calculations were obtained for the frames of 
Figure 4.2.1, where the joints were assumed to be rigid. The first 
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results are for the case when the frames buckled as planar structures. 

( 1 , 2 ) 

These results were computed to compare with existing or known results 
as well as to provide check points for the three-dimensional results. 

The comparison is shown in Table 4.2.2. For these buckling loads, the 
associated sidesway buckling mode shapes are indicated by a dotted 
line in Figure 4.2.1. 

The frame of Figure 4.2.2 together with the indicated loading 
is used as the typical example of the three-dimensional finite element 
buckling theory. For various values of ALPHA (a), buckling loads were 
calculated. Table 4.2.3- Calculations were made for the case where 
the columns and beams were approximated by both one and two elements. 
Using these results, the graph of P CR versus ALPHA, Figure 4.2.3, 
was obtained. Curve I shows results for a one element approximation 
where the columns were not allowed to twist (column GJ = 00 ) . Curves 
II and III are results for one and two element approximations to the 
three-dimensional structure where the columns were allowed to twist. 

The eigenvector values (Curve III), using the procedure described 

in Chapter H are given in Table 4.2.4. With these values, the mode 
shapes for a = 0.0, 0.5, and 1.0 are shown in Figures 4.2.4 (a), (b), 

(c), (d), (e), and (f). These mode shapes are typical examples of 
those obtained for the problems. The mode shapes have been normalized 
so that the u displacements of node 1 are unity. Also, the following 
disp3.acements axe assumed to be equal (inextensional behavior): 

u 2 = U 1B = u l i u 4 = U 3B = u 3 ' Y3 “ W 3C = W 1 i and w 4 = W 2C = V 2 ’ 
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P CE. <“) 

Alpha 

a 

I 

without torsion 
one element 

II 

with torsion 
one element 

III 

with torsion 
two elements 

0.0 

166900 

138980 

138140 

0.1 

162600 

|— — — 

137840 

137040 

0.2 

157800 

136560 

135760 

0.3 

152600 

135060 

134300 

0.4 

147000 

133300 

132560 

0.5 

141200 

131180 

130480 

0.6 

135100 

128580 

127900 

0.7 

129000 

125360 

124740 

0.8 

123000 

121420 

120820 

0.9 

117000 

116720 

116160 

1.0 

-LjJOOO 

r 

111380 

110840 


T ABLE 4.2.3.- buckling loads for three-dimensional, one-bay frame 






























































































































































4.3 Semirigid Joints 


The two-dimensional rigid frame of Figure 4.2.1 (b) is used 
as a semirigid frame as shown, in Figure 4.3 .1. All of the stiffness 
values (R) are the same for each joint. For joints 1 and 3> the 
semirigid joint procedure of Chapter II sis applied to Figure 2.2.1 (a) 
is used. This same procedure, but using Figure 2.2.1 (b) as the 
guide, is applied to Joint 2. The comparison of the buckling loads of 
Figure 4.3.1 with known values based on slope deflection procedures 
with, fictitious members as springs^ ^ is shown in Table 4.3.1* The 
two results are plotted in Figure 4.3-2. Also shown on the plot is 

the range of practical moment connection stiffness obtained by 

(5) 


Rathbun in experimental studies of joints. 


Rathbun's tests 


8 in- lb 

obtained joint stiffnesses that ranged from .001 x 10° — * ■ to 

rad. 

4.5 x 10® for the standard AISC pinned connections and 

rad. 

8 in Th „ ft in-lb. 

from 1.9 x 10 - i a to 180 x 30 for moment connections. 

rad. rad. 

The semirigid, three-dimensional frame. Figure 4.3*3 is 
developed in a similar manner as the previous frame. The results 
from Figure 2.2.2 (a) are us<- o obtain Figure 4.3.3 to 9 Z for 

members 1-2 and 3 - 4, 9 X for members 1-3 and 2-4; and "k° 

9 ). For most connections, the relative torsional slippage at the 

v 


joint will be very slight. 


(*> 


Therefore, the torsional stiffness 


of each member, , was assumed to be constant over the member 

(beam or column). Therefore, torsional semirigid joints are not 
considered in this thesis. The buckling loads for the various values 
of Rj_ and R2 are computed and shown in Table 4.3.2. Using 
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representative values of these eigenvalues, the graphs of versus 
P (Figs. 4.3.4, 4 . 3 . 5 , and 4.3.6) were obtained. This is to 
show the upper and lower limits of each type of connection for a = 0.0, 
0.5, and 1.0. Figure 4 . 3.7 shows the effect of a on the buckling 
load for various values of the Joint stiffness. Also shown is the 
result for the joint stiffness corresponding to the weakest 
stiffness obtained by Rathbun in studies of typical moment connections. 

Figure 4 . 3.8 is a sample plot of load versus determinant 
value at a = 0 for the three-dimensional, one-bay frame. This graph 
is presented to illustrate the procedure used in this paper and 
described in Section 4.1. 


Figure 4.3.1.- Semirigid, two-dimensional* two-bay frame. 


TABLE 4.3. 
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9 rad 
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1.- BUCKLING LOADS FOR SEMIRIGID, TWO-DIMENSIONAL, TWO 
BAT FRAME 
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TABLE 4.3,2.- BUCKLING LOADS FOR SEMIRIGID, THREE-DIMENSIONAL 

' ONE-BAY FRAME 





































































































































Figure 4. 3 . V.- Effect of semirigid Joints on buckling loads of 
three-dimensional, one-bay frame (a = 0.0). 


Joint stiffneafl ft , 10 inch-pounds/radian 












Range of practical moment connection (ref. 5) 



Effect of load parameter and semirigid joints on buckling load* for 
three-dimensional, one -bay frame. 




CHAPTER V 


DISCUSSION OF RESULTS 

The buckling load given in Chapter IV for the two-dimensional 
one-bay frame is the same as for the case when a = 1.0 for the 
three-dimensional problem. The difference from the known literature 
value was only 1.01 per cent. Curve I of Figure 4.2.3 neglects tor- 
sional deformations in the columns and was found to be approximately a 
linear relationship between P and a. Curve II considers torsion and 
shows that the influence of torsion can reduce the buckling load 
by as much as 27,920 pounds or 20 per cent (a = 0.0). Therefore, 
by not allowing the columns to twist, it would be possible to have 
a 20 per cent nonconservative estimate of the buckling load for the 
frame. As expected, the curves (I and II) became the same when 
a = 1.0. The two-element curve (IH) is shown to compliment the 
study of convergence in the same manner as the cantilever beam. 

The two-dimensional, two-bay rigid frame was found to be the 
upper limiting case (R -»Q0 ) for the semirigid one. The semirigid 
buckling loads agree very favorably with the known values. The 
three-dimensional rigid frame corresponds to the upper limiting 
case (a = 1.0 and ~*(D) for the semirigid example. Also when 
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5 

a = 1.0 the same eigenvalue was obtained for Rp = 1 x 10 to 00 . 

This is due to the fact that G = 0 and stiffness Rp is not 

y c 

necessary. Also, as expected, the semirigid buckling load curve for 
a = 1.0 is slightly higher than the semirigid, two-bay frame. 

To develop complete rel:‘ :.c ility of the results presented, the 
eignevalues were checked by the following methods: (l) carefully 

plotting P versus determinant value (Pig. 4*3*®) to eliminate the 
possibility of missing the first buckling load; (2) obtaining 
eigenvectors to assure agreement with the assumed ones; (3) performing 
the check points (upper and lower limits) previously mentioned; 

(4) developing the study of convergence (chapter III) to have a feel 
for the expected accuracy; and (5) comparing semirigid joint results 
with the known values. 

Figures 4.3.2, 4.3.4, 4.3-5, 4-3.6, and 4.3*7 indicate that 
standard moment connections in both two-dimensions and three- 
dimensions give a close approximation to rigid joints. However, 
buckling loads obtained by assuming pinned connections instead of 
the more realistic semirigid joints considerably underestimates the 
buckling loads of the frames studied. 



CHAPTER VI 


CONCLUSIONS 

A general theory for buckling of three-dimensional space frames 
taking into account Joint flexibility has been studied. This 
procedure, while directed toward buckling problems, is also applicable 
for determining stresses and deflections of general finite element 
structures. This method is ideally suited to high-speed digital 
computation by use of a completely automated computer program. 

It has been shown from the convergence study that in all cases 
the buckling loads will converge to the correct answer by increasing 
the number of finite elements. If a high degree of accuracy 
for calculating the eigenvalue is not needed (greater than about 
1 per cent), one element is sufficient to represent the sidesway 
buckling behavior. By comparing the pinned column to the symmetrical 
buckling case, two elements are adequate for computing the buckling 
loads. Results on semirigidity of Joints show that buckling loads 
will be conservative if assumed to be pinned and slightly 
unconservative if assumed to be rigid. 
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APPENDIX 

SAMPLE COMPUTER PROGRAM SHOWING THRffi-DIMHISIONAL SHOKIGID BUCKLING LOAD 

P*F6016. LRC computer complex 

JOB* 1 ,2000, 060000* A 1 774 * ,2, JAMES W RAMSEY ,GMF407, MS 254 , 


OOOOC 4 


PROGRAM 3UCKLG f INPUT, OUTPUT « T APES* 1 NPUT « TAPE6=0UTPUT » 

DIMENSION A 1 3 2, 32 1 ,IM(4)«LG(2), LAMB I 21 • IPlVOTt 321, IR(5>,JNf 21 

000003 


REAL IM,LG,IP*LR»LAMB,JM 

000003 


WRITE 16* 1 1 

000007 

1 

FORMAT 1 1 X* 53H PLfASF SFNO RESULTS TO JIM RAMSEY, MS254, PHONE 4826 


1 

i///lX,36H SIDESWAY BUCKLING LOAD VERSUS ALPHAI 

000007 


READ (5,7) (IM(I), {=1, 4>,(LGf 11,1=1,2) 

000027 

? 

FORMAT ( 6F8 • 1 ) 

000027 


DO 41=1,3 

000031 


GO TO (5,6,7) ,1 

000037 

5 

ALPHA=0.0 

000040 


GO TO 100 

000041 

6 

ALPHA=0.5 

000043 


GO TO 100 

000043 

7 

ALPHA=1 .0 

000045 

100 

00 4K=1 ,10 

000047 


GO TO (30,31, 32, 33,34, 35, 36, 37, 38, 39), K 

000064 

30 

R=l.0E*5 

000066 


GO TO 14 

000066 

31 

R=1.0E*7 

000070 


GO TO 14 

000070 

3? 

R = 3.0F-*-7 

000072 


GO TO 14 

00007? 

33 

R = 8 . OF + 7 

000074 


GO TO 14 

000074 

34 

R = 1 *4 E + 8 

000076 


GO TO 14 

000076 

35 

R=2.1E*3 

000100 


GO TO 14 

000100 

36 

R =3# OF + 8 

000102 


GO TO 14 

000102 

37 

R = 5 * OF + 8 

000104 


GO TO 14 

000104 

38 

8= 1 • 0 ^ + 9 

000106 


GO TO 14 

000106 

30 

R= 1 * 0F+ 10 

000110 

14 

EM= 30* OF + 6 

000112 


GM= 1 1 * 2F + 6 

000113 


JM{ n=o.7?o 

000115 


JM ( 2 ) ^0, 540 

000116 


X = 0*0 

000117 


DO 15N= 1 ,4 



000120 
O00L26 
000127 
000130 
000131 
000140 
00014? 
000146 
000147 
000153 
000155 
000157 
00016 2 
000165 
000170 
000173 
000175 
000176 
000204 
000206 
000214 
000217 
000222 
000223 
000225 
000226 
000227 
000230 
000241 
000243 
000246 
000250 
000251 
000252 
000256 
000257 
000260 
000264 
000266 
000270 
000271 
000274 
000275 


Y=10.0**C6-N1 
00 17H= 1 T 50 
00 16 J = 1 *32 
00 16 L= 1 • 32 
16 A ( J * L ) =0 • 0 

!R(ll=tN(?>/IM(ll 
lR(21=GM*JM(il/(EK*lH( lM 
IR(3)=IM(31/IHC11 
IRI4)=GM*JHI21/{EM*IMI 111 
IR(51-IM(41/IM( 1) 

IR=LG(2 l/LGC 1 1 

C = R*LG(! 1 M 

0=< U0F+9)*LGm/IEM*IMCl U 

7=H/5 • 0 

P=X+Z*Y 

IF <P.E0. 160000. IGO TO 4 
P=-P 

LAMB( I)=P*LG(l»**2/( 30.0*EM*IM<111 
LAH8 ( 2 ) =AL PH A*LAHB (1 1 

AU ,1)=24.0*( 1.0+IR1 3I/LRn72.0*tAMB(ll 
A( 1*4)=6.0*IRI3)/LR 
A( l f 7>=-6*0-3.0*LAMBCl) 

A(till)=ACU4) 

A( 1 t 141 =A( 1*71 
A ( 1 « 1 6 ) =— 4 • O* A < 1 » 4 > 

A(M8I=AU 4 
A C 1 #241 =A ( 1 **t I 

A( 2*2 ) =24.0*1 IRU 1*IRC 3l/LRl+36.0*tLAMe< 1I+LAMBI21 1 
A(2 f 3)=-Al 1*41 

A(2,5)=6.0*IR<1)*3.0*L AMB(l) 

A ( 2 * 9 ) =A( 1 * 16 ) 

A ( 2 * 10 ) = AI 2 * 3 ) 

A(2*17)=A(2*3) 

A( 2*19)=6.0*IR<m3.0*lAN8<2) 

A ( 2 * 2 3 ) =A (2*3) 

A(3*2)=A(2*3) 

A(3?3)=4.0*IR(3I/LR+D 
A(3*9)=A( 1,41 
A(3*10)=A( 1*41/3.0 
A ( 4 * l ) = A ( 1 * 4 ) 

A(4*4)=4.0*IR( 3 1 / L R+ D 
AC4*16»=-A( 1*41 
A ( 4 * 1 8 ) =A (3*101 



000277 

A( 5 » 2) = A( ?»5) 

000300 

A<5,5)=4.0*( IRC ll + LAHBU 1 )*OlR(4)/LR 

000307 

At 5 ,6 > =— C 

000310 

A 1 5 1 1 2 1 IRC4I/LR 

000312 

A t 6 * 5 ) =— C 

000313 

A(6,6)=t>4.0*IR(5)/lP 

000316 

A<6,?0) = 2.0*IP(5)/1_R 

000320 

A(7 ? 1)=A( 1,7) 

000321 

A(7,7)=A<5,5)*4.0*tl.0-IR(l> ) 

000326 

A ( 7, ft ) =-C 

000327 

A( 7,2 1 1 =A( 5,12) 

000331 

A ( 8 , 7 1 =— C 

000332 

A ( ft,ftl=A(6,6) 

000333 

A ( 6 , 1 5 ) = A t 6,20) 

000335 

A( 9,2 ) -A( 2,9 ) 

000336 

A<9,3)=A(1,4) 

000340 

A t 9 ,9 1 =A (2*2) 

000341 

A (9, 10 ) = At 1,4) 

000342 

A(9,12)=A<2,5) 

000344 

A { 9, 17 ) = A( 1*4) 

000345 

A ( 9,23 ) =A ( 1,4) 

000346 

A { 9 , ?5 ) =A (2,19) 

000347 

At 10* 21 =At 2*3) 

000351 

At 10* 3) =A( 3* 1 0) 

000352 

At 1 0*9 1 = A 1 1*4) 

000353 

At 10* 10I=A( 3* 3) 

000355 

A ( 1 1 » i ) =A (1,4) 

000356 

At 1 1 * 1 1 ) =A (4,4) 

000357 

At 1 1, 16)® At 2,3) 

000361 

A( 1 1*24)=A( 3* 10) 

00036? 

A( 12*5 )=A( 5*12) 

000364 

At 12*9) =At 2,5) 

000365 

A( 12, 12) -A (5*5) 

000367 

A ( 12,13) =-C 

000370 

At 13, 12)=-C 

000371 

A ( 1 3* 1 3 ) =At 6 * 6 ) 

000373 

A(i3,26)=A(6,20) 

000374 

A( 14,1 )=A( 1,7) 

000376 

At 14, 14) =At 7, 7) 

000377 

A ( 14,15) =— C 

000400 

At 14,27) =A( 5* 12 ) 

000401 

At 15, 8) =A( 6,20) 

000402 

At 15,14) C 


000403 

000405 

000406 

000410 

000411 

000416 

000417 

000422 

000423 

000425 

00Q4Z6 

000427 

000431 

00043? 

000433 

000435 

000436 

000440 

00044V 

000446 

000450 

000451 

000453 

000454 

000455 

000456 

000460 

000465 

000466 

000467 

000471 

00047? 

000474 

000475 

000477 

000500 

000501 

00050? 

000504 

000505 

0005G7 

000510 

000512 


At 15 • 1 5 )=A( 6,6) 

At 16*1) =A ( 1*16) 

A ( 16*4) — — A ( 1*4) 

At 16* 11 )=A(2* 3) 

At 16* 16 )= A ( 1 tl)*72.0*(LAMBC2 )— LAMB 1 1 ) ) 
A(16*18)=At2*3) 

At 16.21 )— -6.0*3. 0*LAMBt2) 

At 16* 24 ) =A ( 2 * 3 ) 

A 1 16* 27 ) »A C 16.21 ) 

At 17, 2) =A (2.3) 

At 17,9)=A( 1,4) 

At 17, 17 ) =A 13*3) 

At 17, 23) =A( 3. 10) 

A 1 18* 1 ) =A t 1,4) 

A ( 1 8, 4 ) = A ( 4 , 1 8 ) 

At 18,16)^A(2.3) 

At 18, 18 ) = A (4,4) 

A ( 1 9 , 2 ) = A { 2 » 1 9 ) 

At 19,l9)=A(5,5)*4.0*(LANB(2)-LAMBt 1) ) 

A ( 19,20 ) =-C 
At 19*25)-A(5, 12) 

A(20»6)=A(6,20) 

A ( 20. 19 ) =— C 


A ( ?0. 20 ) =A (6*6) 

A(2l*7)=AC5»12) 

A( 21 * 16)=A( 16,21) 

A(21,?1)=A(7,7)+4.0*(LAMB(2)-LAMB( 1) ) 

A ( 2 1 , 2? ) =*C 


A ( 22 1 

,21) 

=-C 

At 

22, 

,22) 

= A ( 6* 6 ) 

At 

22, 

,28) 

=A ( 6 * 20 ) 

At 

23i 

, 2 ) = 

A( 2 , 3 ) 

At 

23, 

, 9) = 

At 1,4) 

At 

23, 

,17) 

= A t 3 , 10 ) 

At 

23, 

,23) 

=A ( 3, 3 ) 

At 

24, 

,1) = 

At 1*4) 

At 

24, 

,11) 

-A ( 3 , 10 ) 

At 

24, 

,16) 

=A ( 2* 3 ) 

At 

24, 

,24) 

—A (4*4) 

At 

25, 

, 9) * 

At 2,19) 

At 

25, 

,19) 

- A ( 5 * 12) 

At 

25, 

,25) 

=A( 19,19) 

At 

25 1 

,26) 

=-c 


000513 

A(26* 13)=A(6*?0) 

000515 

A{26,25)=-C 

000516 

A(26*26)=A<6*6) 

000517 

A( 27* 1 4 ) -A ( 5 * 12 ) 

000521 

A(27*16)=A( 16*21) 

000522 

A(?7,27)=AC 21*21) 

00052^ 

A{ ?7 ♦ 28 ) =-C 

000525 

A(?8*?2)=A(6*20) 

000526 

A(?8*27)=—C 

000527 

A ( 2 8* 28 ) =A (6*6) 

000531 

A ( 3 » 29 ) =-D 

000532 

A ( 4 * 29 ) = — D 

000533 

A ( ?9 t 3 ) =-0 

000534 

A ( 29 * 4 ) =-0 

000535 

AC?9,29)=IR<2 )+2.0*0 

000540 

A( 1 0 * 30 ) ~~D 

000541 

A ( 1 1 *30) --0 

000542 

A(30*10) 0 

000543 

A|30 v in*-0 

000544 

A(30t30)=At29*29) 

000546 

A « 1 7 * 3 l ) =- D 

000547 

AU5,3i)=-0 

000550 

A(31, 17)^0 

000551 

A { 3 1 * 18) =-[) 

00055? 

A(31*31)=AC29*29) 

000553 

A(23* *<?)=-0 

000554 

A(24,32)=-D 

000555 

A ( 32 * 2 3 ) =— 0 

000556 

A ( 32 * 24 ) =“0 

000557 

A( 3? * 32 ) = A ( 29*29) 

000560 

CALL DETEV i A , 32 , DET FR M, 1 PI VOT * 32* I SCALE ) 

000564 

F(OETERM) 15*4*17 

000566 

1 7 Lf-NTINUE 

000570 

r /-XMZ-0.20)*Y 

000577 

*=-X 

000600 

4 Uki TF (6* 18)P*ALPHA,R 

000616 

18 FORMAT (10X*F16.5*F6.2*E16.5) 

000616 

STOP 

000620 

END 


65 


PL EASE S FNO RESULTS TO JIM RAMSFY, MS254, PHONE 4826 


SIOFSrfAY BUCKLING LQAO VERSUS ALPHA 


-4.56400F+04 

0.00 

1.00000E+05 


-9.08800E+04 

0.00 

1.000006+07 


-1 . 15360F+05 

0.00 

3 • OOOOOE+07 


-i .285R0C+05 

0.00 

8.00000E+07 


- 1 . 32 700 E + 05 

0.00 

1.40000E+08 


-1. 34640F+05 

0.00 

2 • 1 0000E+08 


-1 .35840F+05 

0.00 

3.00000E+08 


-1.36980E+05 

0.00 

5.00000E+08 


-1 .37840F+05 

0.00 

l.OOOOOE+09 


- 1 • 3 R640E + C5 

0.00 

1 • OOOQOE+ 1 0 


-4 • 00800F+04 

.50 

i .00000E+05 


-8.42400F+04 

.50 

1 .OOOQQE+07 


-I .08180E+05 

.50 

3 • OOOOOE +07 

■ 

-1.21 120E+05 

.50 

8.00000E+07 


- 1 • 25 160F + 05 

.50 

l .40000E+08 


“ 1 • 27060E + 05 

.50 

2*1 OOOOE+OR 


-1.28220F+05 

.50 

3.00000E+08 


-1 .29340F+05 

.50 

5.00C00F+08 


-1. 30200F+05 

.50 

1.00000E+09 


-1 .30980F+05 

.50 

1 .OOOOOE* 1 0 


— 3 • 22400E+04 

1.00 

1.00000E+05 

4 

-7.06R00E+04 

1.00 

1.00000E+07 


-9* 16000F+04 

1.00 

3.00000E+07 


- 1 • 02 84QF+05 

1.00 

8. OOOOOE+07 

-i 

-1.06300E+05 

1.00 

1.40000E+08 


- 1 • 0794 0F + 05 

1.00 

2.1 OOOOF+08 


-1 . 08940E+05 

1 .00 

3.00000E+08 


-1.09900F+05 

i.oo 

5.00000F+08 

i 

-1.1064QE+05 

1.00 

1 .00000E+09 

4 

-1.11300E+05 

1 .00 

1 • OOCOOE+1 0 

j 

02/19/68 LRC SCOPE 3.0 
14. 26.50. RMF6016. 

6600C- 13 IK 02/12/68 

4 

14.26.50. 

LRC 

COMPUTER COMPLEX 


14.26.50. JOB* 1 *2000.040000 

1 # 

A 1 77 4 * *2* 

1 

14.26.50. JAMES W RAMSFY . 

14.26.51 . RUNC S) 

14.26.58. SETINDF . 
14.26.59.LG0. 

14. 37. 15. STOP 

GMF407* 

MS 254 , 

1 

14.37. 15. CPU 385.822971 

SEC. 


£ 

14. 37.1 5. FPU 006.159488 

SEC. 


:: £ 

RMF6016. PR 1 N T— PP 00365 

LINES 





